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Abstract 

We investigate lowest-lying scalar meson properties predicted from 
QCD Laplace sum rules based upon isovector and isoscalar non-strange 
qq currents. The hadronic content of these sum rules incorporates de- 
viations from the narrow resonance approximation anticipated from 
physical resonance widths. The field theoretical content of these sum 
rules incorporates purely-perturbative QCD contributions to three- 
loop order, the direct single-instanton contribution in the instanton 
liquid model, and leading contributions from QCD-vacuum conden- 
sates. In the isovector channel, the results we obtain are compatible 
with ao(1450) being the lowest-lying qq resonance, and are indicative 
of a non-qq interpretation for ao(980). In the isoscalar channel, the 
results we obtain are compatible with the lowest lying qq resonance 
being fo(980) or a state somewhat lighter than fo(980) whose width 
is less than half of its mass. The dilaton scenario for such a nar- 
rower cr-resonance is discussed in detail, and is found compatible with 
sum rule predictions for the resonance coupling only if the anomalous 
gluon- field portion of dominates the matrix element < d|©^|0 >. 
A linear sigma-model interpretation of the lowest-lying resonance's 
coupling, when compared to the coupling predicted by sum rules, is 
indicative of a renormalization-group invariant light-quark mass be- 
tween 4 and 6 MeV. 

1 Introduction: Status of the Lowest-Lying 
Scalar Resonances 

At present, there is a great deal of confusion concerning both the identity and 
interpretation of the lowest lying 1 = and 1=1 scalar resonances, specif- 
ically the four states denoted in the Particle Data Guide by / (400-1200), 
/o(980),_a (980), and a (1450). The nearness of the / (980) and a (980) 
to a AA threshold has led to a widely held interpretation of these states 
as isopartner AA-molecules [1,2,3], as opposed to light gg-resonances (lin- 
ear combinations of uu and (id states). However, the assumption that these 
states are isopartners and/or A" A"- molecules have both been subject to re- 
cent scrutiny. In particular, Morgan and Pennington [4] have disputed the 
A A interpretation of fo (980). An analysis using the Jiilich model for im 
scattering [5] is compatible with a A A interpretation of / (980), but sees 
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ao(980) as a dynamical threshold effect, as opposed to a true resonance state. 
An even more recent analysis of OPAL data [6] supports the consistency of 
a qq interpretation of / (980). 

Theory is similarly ambivalent regarding the / (980) and a (980) scalar 
resonance states. A QCD sum rule analysis [7] based upon correlation- 
function currents chosen to project out KK- molecule states concludes that 
/o(1500) and / (1710) are better candidates for such multiquark states than 
either / (980) or a (980). Moreover, a very recent coupled channel analysis 
of 7T7T scattering [8] suggests that the /o(980) state may really be two distinct 
S-matrix poles (see also [4]), one corresponding to a KK molecule and the 
other perhaps corresponding to a light qq state. 

Indeed, the determination of the lowest lying qq I = and 1 = 1 states is 
of genuine value as a test of our present understanding of QCD, particularly 
its nonperturbative content. Such lowest lying states, when first compared 
with QCD via sum rule methods [9], were necessarily found to be degen- 
erate, as purely-perturbative and QCD-vacuum condensate contributions to 
scalar-current correlation functions cannot distinguish between 1 = and I 
= 1 channels. Of course, this result was first seen to account for the degener- 
acy of / (980) and a (980) as lowest-lying qq isopartners [9,10]. As has been 
emphasized repeatedly over the last twenty years, however, both scalar and 
pseudoscalar channels exhibit significant sensitivity not only to nonpertur- 
bative field-theoretical effects with infinite correlation length (QCD-vacuum 
condensates), but also to instanton effects, the nonperturbative content of the 
QCD vacuum characterized by finite correlation lengths [11,12,13,14,15,16]. 
The instanton component of the QCD vacuum is known to distinguish be- 
tween 1 = and 1 = 1 scalar (and pseudoscalar) states [12,13]. Such a distinc- 
tion is, of course, quite evident in the pseudoscalar channel's large n — rj mass 
difference, though an understanding of the pseudoscalar 1 = 1 channel nec- 
essarily must take into account the first pion-excitation state because of the 
near-masslessness (and concomitantly reduced sum-rule contribution) of the 
pion [17,18,19]. Similarly, the existence of instanton solutions in QCD nec- 
essarily imposes the theoretical expectation that a similar split occur between 
1 = and I = 1 qq scalar resonance states, with the 1=0 state substan- 
tially lighter than its I = 1 isopartner [12,13]. In this regard, scalar meson 
spectroscopy is a genuine test of QCD. 

Recent activity [20,21,22,23] in re-analyzing old tttt and 7rN scattering 
data has led to the reinstatement of a lowest-lying 1 = scalar resonance that 
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is distinct from the /o(980), conservatively labelled by the 1996 Particle Data 
Group [1] as fo (400-1200). Perhaps of equal importance, assuming a (980) 
and / (980) really are either KK isopartners or other non-gg exotica, is the 
Crystal Barrel Collaboration's recent confirmation of an cto(1450) 1 = 1 scalar 
resonance state [24]. If this isovector state is identified as the lowest-lying I 
= 1 qq-scalar resonance (lattice simulations have shown glueballs below 1600 
MeV to be unlikely), it is important to determine whether the identification 
of its I = isopartner with any portion of the / (400-1200) mass range 
(particularly the high end [4,25]) is compatible with the instanton-generated 
mass difference anticipated from QCD. 

Moreover, the fo (400- 1200) has been widely interpreted to be the a- 
particle signature of chiral symmetry breaking anticipated from NambuJona- 
Lasinio (NJL) dynamics [26], the linear sigma- model (LcrM) spectrum [27], 
as well as in models for qq scattering in an instanton background [12,13] 
and in one boson exchange models of the nucleon-nucleon potential [28]. 
Such a a-particle, however, does not characterize the nonlinear sigma model 
(NLaM); indeed, the empirical absence of a credible a prior to 1996 has pro- 
vided impetus for the development of NLaM ideas into a chiral perturbation 
theory framework. Clearly a clarification of the properties, or even the exis- 
tence [29] , of a light cr-resonance is required to distinguish between LaM and 
NLcrM alternatives for effective theories of low-energy hadron physics. 

We choose here to distinguish, somewhat arbitrarily, four different alter- 
natives for the f (400-1200) resonance that each have some empirical support: 

1. The resonance is both very light {m a < 500 MeV) and very broad (T a > 
500 MeV), as suggested by DM2 data [30] and by Tornqvist and Roos's 
analysis of nn scattering [23] . It must be recognized, however, that such 
a resonance may be generated dynamically [8] and is not necessarily a 
qq state. 

2. The resonance is a-like in mass (500-700 MeV) but substantially nar- 
rower in width (r CT < m a /2), consistent with parameter ranges ex- 
tracted by Svec [20], S. Ishida et al [21], and Harada, Sannino, and 
Schechter [22]. 

3. The lowest I = qq scalar is to be identified with a qq pole at (or 
perhaps masked by) the / (980) resonance [8], with a narrow (Y a < 150 
MeV) width. 
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4. The lowest I = qq scalar resonance is a very broad structure in the 
nix scattering amplitude characterized by a mass comparable to [4] or 
substantially above [25] that of the f (980). 

There is, of course, some blurring of the boundaries between these alter- 
natives, particularly Alternatives 2 and 3. Svec [31] has recently reported 
a single-pole fit of 7r-N (polarized) scattering data leading to a mass (775 ± 
17 MeV) somewhat larger than NJL-LoM expectations, and a width (147 ± 
33 MeV) substantially narrower than those already quoted in support of Al- 
ternative 2. Moreover, Harada, Sannino, and Schechter have demonstrated 
[32] how a very light, very broad state [Alternative 1] can transmute into 
a heavier, narrower state [Alternative 2] when p exchanges are taken into 
account. From a theoretical point of view, Alternatives 1 and 2 support the 
existence of a cr-particle, though straightforward LcrM expectations would 
favour the mass range of Alternative 2 and the broad width of Alternative 1 
[33,34]. Theoretical arguments for a narrower cr-particle more fully consistent 
with Alternative 2 have been advanced through identifying the o with the 
near-Goldstone particle of dilatation symmetry in the strong coupling limit 
[35,36]. Such a dilaton would be expected to have a width similar to that of 
the p-meson, corresponding to cancellation of an enhancement factor of 9/2 
in a naive calculation of the width [33] against an anticipated suppression 
factor l/d 2 a ~ 1/4, where d a is the anomalous mass dimension in the strong 
coupling limit. [] 

In the present manuscript, we employ QCD Laplace sum- rules as a tech- 
nique particularly well-suited to relate the field-theoretical content of QCD to 
lowest-lying resonance properties [14]. We use sum- rule methodology specif- 
ically to address the following questions: 

Which, if any, of the four alternatives discussed above for the lightest I = 
qq scalar are supported by QCD sum rules? In particular, do QCD sum 
rules rule out either of the alternatives that are consistent with an NJL/LcrM 
cr-like object? If the existence of a a is consistent with QCD sum rules, is 
such a a a broad LcrM object, or a narrower strong-coupling dilaton? 

1 Note that appearing in eq. (13) of [35] should be understood to be 131 MeV, not 93 
MeV, so that the d a — 1 prediction of Y a coincides with the LoM-equivalent prediction in 
[34]. This point has been verified through personal communciation with V. A. Miransky. 
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What is the mass range for the lightest I = 1 qq scalar resonance? In partic- 
ular, can we rule out all but exotic interpretations of the ao(980), and does 
there exist sum-rule support for the recently confirmed a (1450) being the 
lowest-lying qq object in this channel? 

In the section that follows, we present the sum rule methodology nec- 
essary to address these questions. Specifically, we show how nonzero res- 
onance widths can be incorporated into the hadronic contribution to QCD 
Laplace sum rules, which are argued to be particularly appropriate for study- 
ing lowest-lying resonance properties. We also demonstrate explicitly how a 
lowest-lying resonance's nonzero width elevates a sum rule determination of 
that resonance's mass. 

In Section 3 we present the field-theoretical content of appropriate scalar 
current sum rules. We discuss the sum rule contribution arising from the 3- 
loop order purely-perturbative QCD contributions to the scalar current cor- 
relation function. We also present nonperturbative sum rule contributions 
arising from QCD-vacuum condensates and direct single-instanton contribu- 
tions to the 1 = and 1 = 1 scalar current correlation functions. 

In Section 4, we utilize the results of the preceding two sections to ob- 
tain a sum-rule determination of the masses of lowest-lying scalar resonances. 
Stability curves are generated leading to estimates of such masses for a given 
choice of width and the continuum threshold above which perturbative QCD 
and hadronic physics are assumed to coincide. Detailed comparison is made 
with earlier sum-rule generated stability curves [9] , showing how the separate 
incorporation of renormalization-group improvement, 3-loop perturbative ef- 
fects, nonzero widths, and the contribution of instantons individually affect 
such curves. 

In Section 5, we examine the isoscalar channel in further detail by ob- 
taining values of the mass, width, continuum threshold, and coupling of the 
lowest-lying qq resonance via a weighted least-squares fit to the overall Borel- 
parameter dependence of the first Laplace sum rule. A relationship between 
the anticipated resonance coupling and a phenomenologically estimable ma- 
trix element is developed in Appendix A. In Section 5 this relationship is 
utilized to obtain an estimate of the light quark mass. This relationship is 
also used to assess the sum-rule consistency of a dilaton interpretation for 
the lowest-lying resonance. 

Finally in Section 6 we present our conclusions concerning the questions 
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we have raised above. We assess the compatibility of sum rule predictions 
for the lowest-lying non-exotic 1=1 resonance with ao(980) and ao(1450). 
We examine in detail the four alternatives presented above for interpreting 
the 1 = scalar resonance spectrum and argue that Alternatives 1 and 4 
appear to be unsupported by a sum-rule based analysis of the lowest-lying 
qq state. We also discuss how our conclusions are affected by possible sum 
rule contamination from higher resonances. 
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2 Sum-Rule Methodology and Lowest-Lying 
(/(/-Scalar Resonances of Nonzero Width 

In the narrow resonance approximation, subcontinuum resonance contribu- 
tions to the light-quark scalar-current correlation function, 2 

n(p 2 ) = i J d 4 xe ip - x < 0\Tj(x)j(0)\0 >, (1) 

j(x) = [u(x)u(x) ± d(x)d(x)]/2, (2) 
are proportional to a sum of delta functions: 



(M(«)) ( 



Im 



= E 



- ~9r 

f (s — m 2 ) + im r r r 
(s — ml) 2 + m 2 r 2 



— ► J2 n 9r$(s-m 2 r ). (3) 

r r -> o r 

The coupling coefficient g r is proportional to m 2 . However, the constant 
of proportionality is expected to be much larger for qq resonances [i.e. res- 
onances that couple directly to the field-theoretical operators in the scalar 
current (2)] than for exotic resonances [37]. It is for precisely this reason that 
sum-rule searches for non-gg scalar resonance states, such as KK molecules 
[7] or glueballs [37,38] , utilize correlation functions based on appropriate KK- 
or gluonic-currents that couple directly to such hadronic exotica. 

Laplace sum rules Rkij) for the scalar current correlation function are 
particularly sensitive to the lowest-lying qq resonance of a given isostructure: 



1 r°° 

R k (r) = - / s k Im[U(s)]e- ST ds 

TV J0 

2 We have normalized 1=0 (+) and I = 1 (-) scalar currents in (2) so as to facilitate 
comparison with ref . [9] . 
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J29rm 2 r k e-^e(s -m 2 r ) 



r 



+ - 



71 Jso 



/ s k Im[(U(s)) pert ]e- ST ds 



(4) 



The summation over resonances in (4) clearly follows from the final line of 
(3). The remaining integral in (4) reflects the anticipated duality [39] be- 
tween purely-perturbative QCD and hadronic physics above some appropri- 
ately chosen continuum threshold s > s . As is evident from (4), higher-mass 
resonances are either absorbed into the continuum [m 2 > s ), or if subcontin- 
uum (m 2 < so), are exponentially suppressed relative to low- mass resonances. 
Consequently, Laplace sum rules are well-suited for determining properties of 
the lowest-lying resonance in a given channel. The subcontinuum resonance 
contribution Rk(r, s ) to the k th Laplace sum rule, defined as 



where mg is the mass of the lowest-lying resonance in a given channel. QCD 
sum-rule methodology for a given resonance channel generally involves ob- 
taining an estimate of the mass rri£ via minimization of the field-theoretical 
content of the left-hand side of (6) with respect to the Borel parameter r 
[14,40]. In practice, the ratio utilized is Ri/R , so as to avoid methodolog- 
ically inconvenient enhancement of continuum and higher-mass resonance 
contributions, as well as heightened dependence on (unknown) higher di- 
mensional condensates, through the factors m 2k appearing in (4). 

The field theoretical content of the first two Laplace sum rules Ro,i(t, s o) 
will be discussed in the section that follows. However, it is important to 
recognize that (6) requires significant modification if the lowest-lying reso- 
nance is broad. If one does not invoke the narrow resonance approximation in 
the final line of (3), but instead assumes a Breit-Wigner (or modified Breit- 
Wigner [1]) shape, one finds that 



I roo 

Rk(r, s ) = R k (r) - - s k Im[(U(s)) pert ]e~ ST ds, 



(5) 



is clearly seen from (4) to satisfy the inequality 




(6) 
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1 m r< s o ,00 a m T 

R t (r,s ) = - E { (s _ m y + m . r ^ * 

= 5]},it' l K 1 r„ r ]mf t -"'e( s ,->na (7) 

r 

with the weighting functions demonstrably related to the narrow reso- 
nance limit (4) via 

\vxLW k (m r ,T r ,T) = 1. (8) 

The net effect of such weighting functions on the lowest-lying resonance con- 
tribution to (6) is to replace m 2 with mfWk+i/Wk. If the lowest-lying reso- 
nance is the dominant subcontinuum resonance in a a given channel, then mn 
can be extracted from the lowest-lying (£) resonance contribution to R\/Rq 
as follows [19,41]: 

W (m e ,T e ,T) ( fli(r,s ) \ 2 

For a given choice of 1^ and s , one can use field-theoretical expressions 
for Ro ; i(t,s ), including both perturbative QCD and nonperturbative QCD 
contributions of infinite and finite correlation lengths [Section 3], to obtain 
from (9) a self-consistent minimizing value of me. This procedure constitutes 
the methodological foundation for the results we obtain in Section 4. 

The weighting functions Wo t i can be derived from a Breit-Wigner reso- 
nance shape by expressing that shape as a Riemann sum of unit-area pulses 
P mr centred at s = m 2 [19]: 

Pm[s - r| = 5SF |e(s " Ml + MT) ~ e(s " M ' ~ Mr)l ' (10) 



n. — >no n *• ' \/ a f 



(s - M 2 ) 2 + M 2 T 2 n ^ V J - / 



ln-j + f r 



j-f 



(11) 



where / is any arbitrarily chosen constant between and 1. If one approxi- 
mates the resonance shape via (11) by truncating n to some finite number of 
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pulses, then the approximation (unlike the n — > oo limit) becomes sensitive 
to the choice of /. The value for / may be chosen to ensure that the area 
under the truncated sum is equal to the area under the "true" resonance 
shape 



f 

j — ( 



1/1 ds ;;. (12) 



/-oo (s - M 2 ) 2 + M 2 r 2 
In an n = 4 approximation, for example, one obtains an area of n by choosing 
/ = 0.70. One finds for this four-pulse approximation [Fig. 1] that [19] 



W k [M,T,r] = 0.5589A fc (m, 3.5119 T,r) 
+ 0.2294A fc (M, 1.4412 T,r) 
+ 0.1368A fc (M, 0.8597 T,r) 
+ 0.0733A fc (M, 0.4606 T,r), (13) 

where 

Afe(M, r, r)M 2k e~ = / P M {s, T)s k e- ST ds. (14) 

J —oo 

In particular, 

a „ \ sinh(MTr) 

MM,T,r)= > , (15) 



and 



Al (M,r,T) = sinh(MTT) 



MTt 



1 

1 + 



M 2 r 



cosh(MTT) 



(16) 



As T — > 0, A k — > 1. Nevertheless, it is easy to show for small values of T 
that [41] 

£j££g = 1 + rV 3 + (r-) (-) 

in which case we find from (13) that W /Wi > 1. 

It is evident from a comparison of (6) and (9) that the introduction of 
a small nonzero width necessarily increases a sum rule driven estimate of 
me [41]. Specifically, if we regard (9) as a constraint implicitly defining a 
function m^(r, T^, so), the result Wq/W\ > 1 for nonzero T implies that 

m £ (r, r>, s ) > ^( r > °) s o)- (18) 
Such behaviour is evident from the analysis carried out in Section 4. 
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3 Field Theoretical Contributions to Scalar- 



Current Sum Rules 



In this section, we seek to identify purely-perturbative and nonperturbative 
QCD contributions to Rq(t) and i?i(r), as defined by (4). The three- loop 
expression for the scalar-current [as defined by (2)] correlation function (1) 
in perturbative QCD (n/ = 3; s = p 2 = —Q 2 ) is given by [42,43] 



If-peri (s) 



3 <-> fa fe 



16tt 2 



l+(^ 

7T 



+ 



a. 



7Y 



95 

45.846 in 

6 



A 17 n 2 
2 ) + 12 



// 



(19) 



The imaginary part of (19), needed for the integrand of (4), is defined con- 
sistent with dispersion-relation conventions: 

2iImU(s) = U{s + ie) -U(s-ie), (20) 
[e.g. Im(£n(—s/fi 2 )) = —n], in which case one finds from (19) that 



Im (U pert (s)) = 



+ 



3s 
16tt 


Mt) 


(t) 


2 


31.864- 










(21) 



We choose to set the renormalization scale via the Borel parameter by setting 
fi 2 = 1/r. Upon substituting (21) into (4), we obtain from (5) the following 
purely-perturbative contribution to Rq(t, sq): 



[R (t, s )] 



pert 



+ 



!6lT 2 T 2 

a 



1 - (1 + s r)e- s " T 
95 fa 



1 + '?)t 



7T 



2 + 

. 3 \ TT 

2 17 rso 



PS T 

I w £n(w)e~ w dw 
Jo 







w [£n(w)} e w dw 



— ) 31.8(3 J 

7T 



(22) 
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As is evident from (4) and (5), the subsequent sum rule Ri(t, sq) can be 
obtained from (22) by explicit differentiation with respect to the Borel pa- 
rameter T. 

Q CD- vacuum condensate contributions to the scalar-current correlation 
function given by (1) and (2) are found via operator product methods to be 
[9,14,44]: 



2 * 

mz 



cond \2Q 2 Q 4 

( 1 7m^ m 2 (Q 2 \\ , 

+ (l6^ + 24^-#(^)) <a ' G * 



2Q 4 2Q 6 
( 27m 2 m 2 ( Q 2 \ , 



+ 0(mf)-(^ + 0(mf)\<a,(q q r>, (23) 

where m q = m u = m^ in the SU(2)-flavour symmetry limit. The final term in 
(23) is obtained via the assumed vacuum saturation of several contributing 
dimension-6 operators [9]. 

Standard dispersion-relationship arguments link the expression (4) for 
Ro(t) with dU/dQ 2 : 

_dn_i r ImU(s) 

dQ 2 ~7ik (s + Q 2 ) 2 1 j 

Equation (24) may be expressed as a Laplace transform (with respect to the 
variable Q 2 ) of the function Rq(t). Since 

-q^ = 1 (re- ST )e- Q T dr ee £(re— ), (25) 



(s + Q 2 ) 
we see from (24) and (4) that 



(26) 
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We note from the Laplace transform definition in (25) that C[t] = 1/ Q 4 , C[l] = 
l/Q 2 , and we find through substitution of (23) into the left-hand side of (26) 
that the leading-order condensate contribution to R is given by [9] 



[Mr)] 



cond 



3 

2 

88tt 



= - < m q qq > + — 
2 167T 



< a s G > 



+ 



27 
0(m q ) 



< a s (qq) 2 > r 



(27) 



We do not include < m q qq > in the order m q terms, as this condensate's 
magnitude (= — is independent of the quark mass [14,45]. 

As discussed in Section 1, both scalar and pseudoscalar correlation func- 
tions are also sensitive to nonperturbative contributions of finite correlation 
length, corresponding to the instanton structure of the QCD vacuum. The 
direct single-instanton contribution to the 1 = 1 pseudoscalar sum rule Rq 
[16,46], 



[Mr)} 



3p 2 



inst. 



1Qtt 2 t 3 



(28) 



is constructed via (4) from a correlation function based on the pseudoscalar 
current 



fix) 



u(x)j 5 u(x) — d(x)~f 5 d(x) /2. 



This direct single-instanton correlator 



(n(p 



2^Y 

inst 



d 4 xe^<0\T[f(x)f(0)] mst \0> 



(29) 



(30) 



can be related to the direct single-instanton contribution to the 1 = 1 pseu- 
doscalar channel of the quark-antiquark scattering amplitude [12] simply 
by tying the external fermion lines together [Fig. 2] to form the vacuum- 
polarization loop corresponding to (30). Consequently, useful properties 
of the single-instanton contribution to qq scattering are also applicable to 
(30). Specifically, the 1 = 1 pseudoscalar qq scattering amplitude in a single- 
instanton background is equivalent to the 1 = scalar qq amplitude [12,13], 
the result of compensating sign changes occurring within the amplitude when 
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«75 ~~ ► 1) and when I = 1 — > / = [12]. These features allow the expression 
(28) to be identified with the instanton contribution to the 1 = scalar- 
channel sum rule, as well as with the negative of the corresponding sum rule 
for the 1 = 1 channel: 

[Ro(r)]- nst = [Ro(r)]l: s ° t = - [Ro(r)]l=] . (31) 

To summarize, the aggregate field-theoretical contribution from QCD to 
the leading Laplace sum-rule (5) is 

[^(r.ao)]^' 1 = [Ro(r,s )] pert 

+ [^o(r)] cond +[ J R (r)]2 1 , (32) 

where the terms on the right-hand side of (32) are respectively given by (22), 
(27), and (31) [via (28)]. Nonperturbative order parameters, specifically the 
condensates appearing in (27) and the instanton-size p in (28), are known 
from other theoretical and phenomenological studies - the values we quote 
as standard in Section 4 are consistent with those of refs. 9, 14, and 16. 

As noted earlier on, the next-to-leading sum-rule Ri(t,s ) can be ex- 
tracted via differentiation with respect to the explicit r dependence of (32): 

d 

Ri(t, s ) = -—Rq(t, s ). (33) 

However, the expressions Ro(t, s ) and Ri(r, s ) obtained from (32) and (33) 
require renormalizat ion-group (RG) improvement to be useful - otherwise 
results obtained via (9) are unnaturally dependent on the specific choice for 
a s . QCD Laplace sum rules have been shown to satisfy RG equations 
with respect to the Borel scale parameter r [47]. Consequently, the sum- 
rules i? ,i, once obtained from (32) and (33), can safely be RG- improved by 
replacing a s with the running coupling-constant a s (r -1 / 2 ). In the Section 
that follows, factors of a s in i?o,i will be understood to correspond to the 
PDG [1] 3-flavour expression for a s (r -1 / 2 ) to 3-loop order. 



3 Such is the case in ref. 9, where a s is chosen to be 0.6. 
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4 Sum Rule Analysis of qq Scalar Resonances 



For a given choice of 1^, so, and r, one can obtain a self-consistent value of the 
lowest-lying resonance mass mi by solving (9), as noted in Section 2, thereby 
defining a function m^(r, Tg,so) implicitly. Such a procedure is meaningful 
provided only the lowest-lying qq resonance dominates the sum-rule. This 
assumption is a reasonable one provided other subcontinuum resonances in 
the same channel are either much heavier, and thus exponentially suppressed 
(4), or are non-gg exotica, as has already been noted. A clear signal of this 
assumption's validity is the absence of a (980), if exotic, from the sum rule 
generated from 1=1 gg-currents (2), even though a (980) is the lowest lying 
1 = 1 scalar resonance. One cannot explain such an absence if ao(980) is 
qq; however its observed absence (see below) is consistent with it being both 
exotic Q and sum-rule decoupled. 

For a given choice of Yi and so, we shall examine the r-dependence of 
m^(r, T^, so) obtained in a given channel from (9) via the field-theoretical 
expressions for -Ro,i( r , s o) developed in Section 3. For such externally imposed 
values of and Sq, we shall identify the true value of mi with the minimum 
value of m^(r) obtained over an appropriate range of the Borel parameter r = 
1/M 2 . Since M = r- 1 ' 2 is itself the renormalization-scale (fi) for the field- 
theoretical content of Ro,i(t, so), a sum-rule calculation cannot be meaningful 
unless we restrict this mass scale to be not only well-above Aqcd, but also 
to be bounded from above by the continuum threshold: M 2 < sq- This 
criterion is necessary to ensure that the continuum contribution to (22), 
which increases as Sqt = Sq/M 2 becomes small, is not overly large compared 
to the remainder of the purely perturbative contribution, i.e. the Sq — > oo 
limit of (22). A more realistic criterion, in which continuum effects are less 
than 30% of the purely perturbative contribution [14], would lead to an even 
tighter upper bound on M. Moreover, the true value of rri£ should not just 
be the global minimum of m^l/M 2 , Tg, s ) with respect to M chosen in this 
subcontinuum range; it should also exhibit insensitivity to small changes in 
the Borel parameter M - a local minimum. This criterion of flatness, as 
well as the requirement for a sensible range of the Borel/renormalization 
scale M, is quite standard in sum rule applications [40]. We also employ 

4 i.e.., a K X-molecule [2,3], a dynamical-threshold effect [5], or possibly an ss state 
with Zweig-Okubo-suppressed coupling to the nonstrange current (2). 
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standard values [9,14,16] for QCD's nonperturbative order parameters in our 
analysis: < m q qq >= -/ 2 m 2 /4 (/ ff = 131MeV), < a s G 2 > = 0.045 GeV 2 ,< 
a s {qqf >= 0.00018 GeV 6 , p = (600 MeV)" 1 . Factors of a s that are not 
absorbed in (approximately-) RG-invariant QCD-vacuum condensates are 
replaced with 3-loop rif = 3 running (Aq C d = 150 MeV) coupling constants 
a s (M)[M = r -1 / 2 ], as discussed at the end of Section 3. 

4.1 1 = Scalar Channel: Narrow Resonance Limit 

In Figure 3, the Borel-scale dependence of the lowest-lying 1 = scalar res- 
onance mass (denoted as M a ) is displayed for a number of choices of so, 
assuming zero resonance width. The curves displayed are restricted to values 
of M 2 < so- As in ref. [9], we see that the curves M a (M) increase with 
the choice for sq. We also see that a subcontinuum local minimum does not 
develop until we consider values of sq larger than 1.6 GeV 2 . Figure 4 displays 
explicitly the Sq dependence of the values of M a that are local minima for 
each choice of Sq in Fig. 3. As is evident from Fig. 4, the minimum value 
of M a for a given choice of Sq also increases with sq. Our analysis finds the 
onset of a local minimum occurring when Sq = 1.61 GeV 2 , corresponding to 
M a = 680 MeV; the Figure 4 curve begins with this point. 

It is useful to compare these results to the seminal 1982 sum-rule analy- 
sis of ref. [9], which included the condensate contributions (27), but which 
did not include the then-unknown instanton contribution (28,31) or correct 
higher-order perturbative contributions. f\ Figure 5 demonstrates the cru- 
cial role instanton contributions play in lowering the lowest-lying I = qq 
resonance mass, though at the price of diminishing the broad range of Borel- 
parameter stability observed in [9] at comparable values of Sq. The use of 
RG-improvement also has a significant lowering effect on the lowest-lying I = 
qq resonance mass. In Figure 6, we compare the stability curves obtained 
from the full field-theoretical content of Ro t i(r, s ) at s = 1.55GeV 2 with 
and without RG-improvement of [-R ,i(t, s )]p ert . The upper curve is obtained 
with a s = 0.6, as in [9], while the lower curve utilizes a s (M). Finally, we note 
that instanton and RG- improvement effects are offset somewhat by higher- 
order perturbative contributions. Figure 7 demonstrates how going from two 

5 The analysis considered only 0{a s /ii) 2-loop contributions to (21), with 17/3 erro- 
neously given as 13/3. The analysis also did not incorporate RG-improvement of a s . 
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to three loops can increase the estimate of the lowest-lying 1 = resonance 
mass by as much as 150 MeV, suggesting higher-order effects as a clear source 
of theoretical uncertainty in scalar-channel sum-rule methodology. Since the 
leading 4-loop contribution to Im[U pert (s)} is both large and the same sign 
[43] as the leading 3-loop [O(a^)] contribution to (21), we can anticipate that 
sum-rule determinations of M based upon (21) will likely underestimate M a , 
a point that will prove important in assessing the viability of the broad, very 
light a (Alternative 1) discussed in Section 1. Q 

4.2 1 = Scalar Channel: Nonzero Widths 

Figures 8-12 present stability curves analogous to Fig. 3 that are obtained 
from (9) with input values of Tt between 100 and 500 MeV. The weighting 
functions Wo,i are obtained via the four-pulse approximation leading to (13) 
and (14). For each value of Ti considered, it is possible to construct a (Fig. 
4 analog) plot of the minimizing value of M a for each choice of Sq. These 
plots are presented in Fig. 13, and they clearly demonstrate how sum-rule 
determinations of the lowest-lying mass increase with increasing resonance 
width. 

In comparing Figs. 8-12, we note that the onset of a local minimum 
below the continuum threshold is itself width-dependent. Specifically, the T 
= 500 MeV stability curves of Fig. 12 do not develop a locally flat minimum 
below So until so is larger than 1.8 GeV 2 , corresponding to a minimizing 
value of M„ at onset that is substantially above 1 GeV. Local minima for 
larger values of s Q are seen to lead to even larger values of M a , as is evident 
from the topmost curve of Fig. 13. The values of so and M a corresponding 
to the onset of local minima are tabulated in Table I for < T < 500 MeV, 
corresponding to the initial left-hand points of the six curves displayed in 
Fig. 13. The values of M a in Table I correspond to the lowest sum-rule 
determinations of M a possible for each choice of V that are consistent with 
the methodological constraints delineated at the beginning of this Section. 
The Table entries clearly indicate that a light lowest-lying I = qq scalar 
resonance cannot have a width larger than half of its mass. [] 

6 Our somewhat low value for Aqcd is similarly motivated to bring down the size of 
perturbative contributions, although in practice our results are virtually unaffected by a 
fOO MeV increase in this parameter. 

7 A qualitatively similar conclusion is stated in ref . [38] . 
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This conclusion is softened only slightly if we relax the flatness require- 
ment and insist only on identifying M a with the global minimum over the 
Borel-parameter range (Aq C d) 2 < M 2 < s . As T increases from 100 to 500 
MeV, Figs. 8-12 still show such a global minimum to be increasing with F. 
When T = 100 MeV, such a minimum occurs at about 530 MeV [see the 
So = 1-2 GeV 2 curve of Fig. 8], but as T increases to 500 MeV, this global 
absolute minimum increases to 840 MeV [see the Sq = 1.4 GeV 2 curve of Fig. 
12]- 

These qualitative conclusions concerning the non-viability of a very light, 
very broad qq resonance are not expected to be altered by further refine- 
ments in the treatment of the lowest-lying resonance shape. Once can obtain 
"improved" weighting-factors Wo t i, as opposed to those (13) based upon a 
four-pulse approximation of the Breit-Wigner resonance shape by increasing 
the value of n used to truncate the Riemann sum (11). Figure 14 demon- 
strates that such an increase in the number of pulses used to approximate the 
Breit-Wigner shape does not appreciably alter the local-minimum value M a 
obtained for a given choice of so and T. It should be noted, however, that Mo- 
increases slightly as n increases, suggesting that a more precise modelling of 
the Breit-Wigner resonance shape would only serve to increase the M a values 
of Table I. 

A more fundamental issue is whether the Breit-Wigner shape is appro- 
priate at all for the modelling of broad resonances, as the Breit-Wigner tail 
will extend significantly into Euclidean (s < 0) and continuum (s > s ) re- 
gions if T is sufficiently large. We have relied upon an admittedly crude n=4 
(4-pulse) approximation in order to minimize such unphysical contributions. 
Larger choices of n increase the sensitivity of Wq,i to unphysical regions in s 
by including square pulses of greater width than the largest pulse in Fig. 2. 

It has already been argued that higher-order perturbative contributions 
will, if anything, increase the M a values of Table I. Consequently, the only 
means compatible with QCD sum-rule methodology that is available for ob- 
taining lower sum-rule estimates for M a is to increase the size of the instanton 
contributions [c.f., Fig. 5], possibly via a decrease in the value for p in (28). 
To do so, however, has ramifications for the 1 = 1 scalar resonance channel, 
as discussed below. 
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4.3 1 = 1 Scalar Channel 



We have seen in Section 3 that the sign of the instanton contribution to i?o,i 
in the 1=1 scalar channel is reversed from that the 1 = scalar channel. 
Consequently instanton effects are now seen to increase the scale of sum- 
rule determinations of the lowest-lying contributing 1 = 1 resonance. The 
stability curves presented in Figs. 15 and 16 clearly indicate a much higher 
range of values for masses of such lowest-lying resonances, as well as a need 
to go to much higher values of so in order to attain subcontinuum stability 
with respect to the Borel-parameter M. 

Table II lists values for the lowest-lying contributing 1 = 1 resonance mass 
(Ms) associated with the onset of a subcontinuum local minimum. Since 
local-minima Ms increase with the choice of Sq, the Ms values in Table II, 
like the M a values in Table I, represent the minimum such values obtainable 
via sum-rule methodology outlined in the beginning of this Section. Table II 
clearly indicates a lowest-lying contributing 1 = 1 scalar resonance in excess 
of 1.49 GeV. 

The much larger 1 = 1 values for the continuum threshold in Table II 
are comparable to values for s obtained via sum rule analysis of the 1 = 1 
pseudoscalar channel [18,19]. In both 1=1 channels, sum rule methodology 
would suggest that Borel stability not occur for values of s^ 2 that are less 
than the masses of contributing subcontinuum resonances. The large value 
of Ms evident from Table II necessarily requires values of sq larger than M|, 
consistent with the values of sq actually listed. 

As noted at the very beginning of this Section, the results of Table II 
clearly rule out a contribution from a (980), the lowest-lying scalar resonance 
in the 1=1 channel. Consequently, the coupling of this resonance to the 
gg-current (2) [with a minus sign chosen for 1=1] must be suppressed, as 
would be anticipated form a non-gg interpretation of this state. The results 
of Table II, however, appear compatible with a (1450) being identified with 
the lowest-lying 1=1 qq-scalar resonance, particularly if the width of this 
resonance is less than or of order 200 MeV. 

As a final comment, we note that any attempt to obtain sum-rule support 
for a broad, very light qq state by enhancing the magnitude of the instanton 
contribution (e.g. by decreasing the instanton size p) will necessarily drive 
up the mass of the lowest-lying contributing 1 = 1 state, i.e., the isovector 
qq state. Thus the price of tuning instanton effects to accommodate a light 
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broad isoscalar qq resonance is likely to be an exotic interpretation for both 
a (980) and ao(1450), the two lowest-lying 1=1 scalar resonances. 



21 



5 Global-Fit of Lowest-Lying 1 = Scalar 
Resonance Properties 



In Section 4, we have utilized a specific choice of the Borel parameter corre- 
sponding to minimization of the sum-rule- derived mass. However, Leinweber 
[48] has stressed the value in an overall fit of the Borel-parameter dependence 
of a sum rule's field-theoretical content to the dependence anticipated from 
resonance properties. For example, properties of the first pion-excitation 
state IT have been obtained by fitting the QCD Borel-parameter dependence 
for the lowest LT-sensitive Laplace sum rule to its corresponding hadronic 
content [18,19]. In this section, we will apply this same procedure to the 
scalar 1 = channel. Specifically, we will obtain a x 2 -minimizing weighted 
least-squares fit of [Ro(r, s )] I=0 to the r-dependence anticipated from the 
corresponding resonance contribution (7), assuming the contribution from 
the lowest- lying resonance dominates the 1=0 scalar channel. 

We begin first by modifying the scalar current (2) to include the non- 
strange quark mass in the SU(2) /-invariant limit: 



where m q = (m u + m d )/2. Although the additional factors of the quark 
mass cancel out in an R\/Rq determination of resonance properties, such 
as that in Section 4, these factors have nontrivial consequences when the 
overall r-dependence of Rq itself is being analyzed. The current (34) now 
corresponds to an RG-invariant operator, and upon RG-improvement, the 
sum rules devolving from that current (as in Section 3) will now include the 
additional r-dependence associated with the running quark mass: 




(34) 



[R (t, s )] 



7=0 




3 



{ [l - (1 + s r)e- s " T ] [1 + 17a s (r)/3n] 



167r 2 r 2 



+ 




(35) 
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Q ' M = 9ER hsI^J^H' (36) 



m q (r) 



rn 



[L(r)]V9 
Ut) 



0.1989 - 0.1756 £n[L(r)] 



(37) 



-^n (rA^ D ) . (38) 

The condensate and instanton contributions to (35) are given respectively 
by (27) and (28). The running coupling (36) and mass (37), as well as the 
purely perturbative contribution to (35), are given only to two- loop order 
to accommodate the numerical limitations of our fitting procedure, which 
involves generating resonance parameters via a Monte-Carlo simulation of 
uncertainties (see below). The parameter rh in (37) is the RG- invariant 
nonstrange quark mass. 

Our fit is generated by obtaining values for m a , T a , s , and the resonance 
coupling g a that minimize a least-squares fit of (35) to the lowest-lying res- 
onance contribution anticipated via (7): 

[Ro(r, s )} 1=0 = gvW (m a , r CT , r)e- m ° T . (39) 

In (39), the weighting factor Wo is given by the four-pulse approximation 
expression (13), with A given by (15). Our weighted least-squares fit is 
over the Borel-parameter range 0.4 GeV -2 < r < 2.2 GeV~ 2 [0.67 GeV 
< M < 1.6 GeV]. The region in r is obtained by requiring an uncertainty 
of less than 20% on the theoretical contribution (35) to R , based upon a 
30% continuum and a 50% power law uncertainty. Parameter values for 
perturbative and nonperturbative scale parameters are as given in Section 4. 
As mentioned above, uncertainties in the fitted parameters {m^r^so,^} 
are obtained only for the two loop case via a Monte- Carlo simulation of the 
power law and continuum uncertainties described above, as well as a 15% 
variation in the size of p, and a factor of 2 "vacuum-saturation uncertainty" 
in the value of < a s (gg) 2 >. The parameter values for {m a ,T a , so, g a } we 
obtain correspond to the fit of (35) to (39) that minimizes a x 2 weighted for 
the previously mentioned continuum and power law uncertainties. 

8 Treatment of this quantity as well as the specific form of Monte-Carlo modelled un- 
certainties is discussed at length in ref. [18]. 
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Figure 17 demonstrates the success of these fitted parameters in matching 
theoretical (35) and hadronic (39) expressions for Rq. The results of this fit 
are as follows: 



m a = 0.93±0.12GeV, (40a) 

< T CT < 0.26 GeV, (40b) 

s = 3.20 ± 1.20 GeV 2 , (40c) 

9a = m 2 ■ (0.039 ±0.014 GeV 4 ), (40d) 



where the quoted uncertainties reflect 90% confidence levels. The results 
(40) are clearly indicative of a lowest-lying gg-resonance that is neither very 
broad nor very light. These results are most consistent with the Alternatives 
2 and 3 delineated in Section 1. 

Pertinent to Alternative 3, the results (40) are certainly consistent with 
identifying the lowest-lying 1 = scalar nonstrange qq state with / (980), as 
has been suggested by a very recent OPAL analysis of Z° decays [6]. Indeed, 
the very broad region of Borel-parameter stability characterizing the Sq = 2.4 
GeV 2 curve of Fig. 3 [which yields an / (980) lowest-lying resonance mass] 
corroborates parameter values for m a , T a , and s within the fitted range (40). 

However, the results (40) do not exclude an Alternative 2 qq state some- 
what below the / (980) in mass, consistent with a KK interpretation of 
/o(980). Such an / (980) would be expected to be essentially decoupled 
from sum rules based upon the current (34), as has been discussed in Sec- 
tion 4. This Alternative 2 qq interpretation becomes particularly interesting 
upon examination of the fitted parameter g a , which cancels entirely from 
the R\ /.Ro-based analysis of Section 4. Arguments within a LaM context 
(buttressed by more general PCAC arguments) have been advanced for ob- 
taining the matrix element by which a scalar current couples a physical a to 
the vacuum [49]: 

< a\m q (u(0)u(0) + J(0)d(0)) |0 >= / w m 2 . (41) 
In the above expression, f n is 93 MeV. In Appendix A we show that 

g a = \<a\j s (0)\0>\ 2 . (42) 
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Making use of the 1 = scalar current j s defined via (34), we find from (41) 
that 

9a = fWjA. (43) 
Comparison of (40d) and (43) suggests that 

m = 4.6±J;| MeV, (44) 

a value somewhat on the low side of the expected range [1] for the nonstrange 
current quark mass. It should be noted, however, that a reasonable light- 
quark mass cannot be obtained in a dilaton scenario [35,36] in which 

< <7|6£(0)|0 >= f a ml > /X, (45) 

unless the contribution to the matrix element from the scalar-current com- 
ponent of 9^ is negligible compared to the anomalous gluonic-field contri- 
butions. If one assumes that ~ m q {uu + dd), the quark mass one would 
obtain by replacing (41) with (45) [via comparison of (40d) with (42)] will 
be larger than (44) by a factor of rn^/m^. 

As a final comment, it should be noted that the value for sq in (40c) 
is sufficiently high to exclude / (1370) and / (1500) from the continuum. 
Although such resonances (if subcontinuum) are expected either to be expo- 
nentially suppressed [via (7)], or if non-gg exotica, to be essentially decoupled 
from a sum rule based on qq currents, there is nevertheless reason to be con- 
cerned about the validity of assuming (as we have done in the section) that 
only one resonance contributes to R . To address the possibility of contri- 
butions from more than one resonance, we have examined whether the field 
theoretical content of [Ro(t, s )] /=0 might be better fitted by contributions 
from two subcontinuum resonances by replacing (39) with 

Ro(i~, s )] I=0 = giWoim,!, I\, r) exp(-mir) + g 2 exp(-mlr). (46) 

The second subcontinuum resonance [presumably f (980)] is assumed to 
be narrow: Wo(m2,0,r) = 1. The remaining resonance parameters are de- 
termined via a weighted least-squares minimization of x 2 - We have found 
such minimization to be accompanied by equilibration of m x and m 2 to the 
previously fitted value (40a), indicating that the r-dependence of the field- 
theoretical content of Rq [Fig. 17] is best fit by a single exponential (39), 
rather than by the sum of two distinct exponential contributions (46). 
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This result indicates that the leading 1=0 scalar-current sum rule is dom- 
inated by a single resonance, as has been assumed in the body of this section. 
Such behaviour is to be contrasted with that of the leading 1=1 pseudoscalar- 
current sum rule, for which a fit including two contributing resonances [it and 
7r'] leads to a x 2 an order of magnitude lower than that of a single resonance 
fit [18]. 
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6 Conclusions 



In this paper we have sought to determine which of several empirically jus- 
tifiable interpretations of the scalar mesons are compatible with theoretical 
constraints based upon QCD sum-rule methodology. 

For the isovector channel, we have found in Section 4 that a (980) is 
decoupled entirely from the isovector sum rule based upon the scalar current 
correlation function given by (1) and (2), a result requiring an exotic (non- 
qq) interpretation for this resonance. The Table II values for the mass of 
the lowest-lying resonance in this channel that does couple to the sum rule 
appear compatible with the lowest lying qq resonance being a (1450). 

For the isoscalar channel, we have delineated in the introductory section 
of this paper four Alternatives for the lowest-lying qq scalar resonance. We 
found in Section 4 that the lowest-lying I = qq scalar resonance cannot be 
both very light and very broad [Alternative 1], provided this resonance is the 
sole contributing resonance to this channel's leading QCD sum rules. Thus, 
the only way the results of Section 4 could accommodate such a resonance 
would be if this lowest-lying resonance's contribution is itself masked by 
higher- resonance contamination of (7), the hadronic side of the sum rule. 
The results at the end of Section 5 show this to be unlikely. Unless we 
require the Alternative 1 scenario to specify /o(980) to be an additional light 
qq excitation, as opposed to a KK state, and unless the lowest-lying qq state 
is unnaturally decoupled, such contamination is itself possible only for values 
of So large enough to include more massive resonances in the subcontinuum 
region, such as / (1370) and / (1500). The values of s Q appearing in Table I 
are clearly below this threshold. 

The results of Section 4 show demonstrable support for Alternative 2, a 
narrower and more massive lowest-lying resonance that is still lighter than 
the (presumably exotic) / (980), particularly if the lowest lying resonance 
is the sole contributing resonance to the sum rule. Moreover, the results 
listed in Table I correlate masses for such a resonance with values of s that 
preclude any contamination of the hadronic side of the sum rule by / (1370) 
and / (1500). These results clearly favour masses above 700 MeV and widths 
below 300 MeV, and are suggestive of Svec's most recent single-resonance fits 
to tcN data [31]. 

It should be emphasized that the results of Section 4 do not exclude the 
possibility that / (980) is itself the lightest I = qq scalar resonance. In- 
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deed, this Alternative 3 scenario is strongly supported by the fit performed 
in Section 5, which gives results (40) consistent with / (980) but problemat- 
ical for a dilaton-explanation of Alternative 2. Specifically, the coupling g a 
obtained in (40d) is much smaller than that anticipated for a dilaton, un- 
less the anomalous gluonic component of G>^ dominates the matrix element 
(45). Moreover, sum-rule contamination from higher 1 = scalar resonances, 
though certainly possible for the fitted range (40c) of s , would be expected 
to increase rather than diminish the apparent value of g a in a single reso- 
nance fit, f\ leading to even a smaller true value for g a and an even larger 
discrepancy from the scale (45) anticipated from a dilaton interpretation. 

We therefore conclude that the support QCD sum rules provide to Al- 
ternative 2 [a lowest lying resonance below /o(980) whose width is less than 
half its mass] does not extend to a dilaton interpretation of this state, un- 
less the matrix element (45) is driven almost entirely by O^'s anomalous 
gluon piece. As a final comment, the Alternative 4 scenario, in which the 
lowest-lying qq isoscalar resonance is broad and at least as massive as /o(980), 
clearly contradicts the QCD sum rule results of Section 5. Moreover, even 
if higher-mass resonance contamination of (7) were to occur, such additional 
hadronic contributions would be expected to drive the apparent mass of the 
lowest-lying state above its true value, as has already been noted. Thus the 
fitted result (40a) for the lowest-lying resonance mass in the absence of such 
contamination cannot be reconciled to a true value much in excess of 1 GeV. 



Acknowledgements 

V. Elias and T. G. Steele would like to acknowledge support from the Natural 
Sciences and Engineering Research Council of Canada. V. Elias is grateful 
for detailed discussions with V. A. Miransky and M. D. Scadron, as well as for 
useful information provided by T. Barnes, G. Lafferty and C. Strassburger. 



'Provided the g r in (7) are positive-definite: e.g., (42). 
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Appendix A: The a Contribution to the 1 = 
Scalar Correlator 



To obtain the explicit a contribution to the 1 = scalar-current (j s ) corre- 
lation function (1), we begin by noting that insertion of a complete set of 
intermediate single-particle a states implies that 



<0| T 3s (x)j s (0)\0 > 

- e ( -.)<o W o)(/S^f), W ,o> 

+ (higher-resonance and multiple particle terms) (A.l) 

We will define the matrix element connecting a physical (on-shell) a to 
the vacuum via a scalar current to be 

< 0\j s (x)\a(q) >= Me-'fv^ 1 "-^), (A.2) 

where the Heaviside step function Q(x ) can be expressed as the following 
integral: 

1 poo pixQT 

e(x )=Um— / — dr. (A.3) 

e-*o 2rn J-oo r — i\e\ 

When xq > 0, the contour along the real r axis can be closed by an infinite 
arc in the upper half plane, enclosing the residue at r = i\e\ : Q(xo) = 1. 
When xq < 0, the contour must be closed in the lower half plane, and the 
residue is not enclosed: 6(x ) = 0. If we substitute (A.2) and (A.3) directly 
into (A.l), we obtain 



<0\Tj s (x)j s (0)\Q>* 



-i\M\< 



lim 

e-o (2tt) 4 



/oo rc 
-OO J — i 



dT 



2(r — i\e\)J(p + m 
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+ 



(p+m%-T)xo -iq-x 



2(r - i\e\)J(p + ml 



(A.4) 



where the a subscript on the left-hand side of (A.4) represents the explicit 
contribution of single-cr intermediate states to the matrix element. To eval- 
uate (A.4) further, we make the following change of variable from r to q : 



Qo = +m 2 a -T , dq = -dr. 
We then obtain (d A q = d 3 q dqo ; q • x = go^o — q- x)\ 



(A.5) 



<0| Tj s (x)j s (0)\0 > a 



+ 



lim 



i\M\ 2 



o (2tt) 



d 4 q 



-iq-x 



2(yjq* + m 2 a -q - iM)^ + m 2 



iq-x 



= lim 



+ 



2( v /cf + m 2 a - q - ileD^q 2 + m 2 

i\M\ 
6 (2tt) 



OO p 

a g 



-iq-x 



2^J(p + ml _^J(p + ml 



^Jq 2 + m 2 a + q - i\e\ 



(A.6) 



The final line of (A.6) is obtained from the intermediate line by letting g M — > 
— g M in the second portion of the integral. We find trivially from (A.6) that 
(q 2 = q 2 -q 2 ;e' = 2^ ( ? + ml e); 



<0\Tj s (x)j s {0)\0> a 

i\M\ 2 {<*> e -iq-x 

= hm +t — / d q — — r . (A. 7) 

e '-t) (2tt) 4 7-oo y g 2 -m2+i| e '| V ; 

Upon substitution of (A. 7) into (1), we find that 

(hm)„ = - 'r'^T-M < a - 8 > 
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where we have utilized the explicit definition of the constant M in (A. 2). 
The propagator denominator in (A. 8) can be modified in the standard way 
[1] for incorporating relativistically a Breit-Wigner width T a : 



dV))„ = (A.9) 

pZ _ m z _|_ ima Y a 



Consequently, we see from comparison to (3) that the coupling coefficient g r 
is given by 

g r =\< cr\j s (0)\0 > | 2 (A.10) 
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r(MeV): 100 200 300 400 500 
So(GeV): 1.61 1.61 1.62 1.65 1.70 1.82 
m CT (MeV): 680 687 716 778 884 1087 

Table I: Sigma masses associated with the onset of a subcontinuum local 
minimum in the isoscalar channel. 



r(MeV): 100 200 300 
s (GeV 2 ): 3.24 3.26 3.38 3.62 
M 5 (GeV): 1.49 1.50 1.55 1.63 

Table II: Lowest-lying masses (Ms) associated with the onset of a subcon- 
tinuum local minimum in the isovector channel. 
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Figure Captions: 



Figure 1: An example of the 4 square-pulse approximation to the Breit- 
Wigner resonance shape obtained by truncating eq. (11) to n=4, and 
by choosing f = 0.701 to ensure that the area under the four pulses 
is equivalent to the total area under the Breit-Wigner curve. This 
particular example is for a mass M = 680 MeV and width T = 100 
MeV. 

Figure 2: 

(a) The single-instanton contribution to the quark- ant iquark scattering am- 

plitude (c.f. fig. 5 of ref. [12]). 

(b) The single-instanton contribution to the corresponding current-current 

correlation function. 

Figure 3: Stability curves for determining via (9) the lowest-lying 1 = 
resonance mass. The corresponding width of this resonance is assumed 
to be zero. The stability curves are truncated to exclude the region in 

which the Borel mass-scale M exceeds s l J 2 . 

Figure 4: The s -dependence of the masses M a corresponding to (local) 
minima of the stability curves presented in Fig. 3. 

Figure 5: How instanton contributions affect M a : a comparison of two 
stability curves for the lowest-lying 1 = resonance mass obtained by 
including or not including the direct single-instanton contributions to 
the field-theoretical expressions for _R ,i( r ? s o)- Both curves are ob- 
tained assuming zero resonance width and a continuum threshold so = 
1.55 GeV 2 . 

Figure 6: How renormalization-group improvement affects M a : a compar- 
ison of two stability curves for the lowest-lying 1 = resonance are 
obtained by including or not including the r dependence of the strong 
coupling. The constant coupling curve is obtained by choosing a s = 0.6, 
as in ref. [9]. The (3-loop) running coupling curve is obtained assuming 
Aqcd = 150 MeV. Both curves are obtained assuming zero resonance 
width and a continuum threshold So = 1.55 GeV 2 , and both curves in- 
corporate direct single-instanton contributions to the field-theoretical 
expressions for -Ro,i( r ! s o)- 
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Figure 7: How three- loop perturbative contributions affect M a : a compar- 
ison of two stability curves for the lowest-lying 1 = resonance mass 
obtained by including or not including three-loop order perturbative 
contributions to the field-theoretical expressions for Ro,i{t,sq). As in 
Figs. 5 and 6, s Q is assumed here to be 1.55 GeV 2 , and the resonance 
width T is assumed to be zero. 

Figure 8: Stability curves for determining via (9) the lowest-lying 1 = 
resonance mass, assuming a resonance width of 100 MeV. As in Fig. 
3, the stability curves are truncated to exclude the region in which the 

Borel mass scale M exceeds Sq 2 . 

Figure 9: Stability curves for determining the lowest-lying 1 = resonance 
mass, assuming a resonance width of 200 MeV. 

Figure 10: Stability curves for determining the lowest-lying 1 = resonance 
mass, assuming a resonance width of 300 MeV. 

Figure 11: Stability curves for determining the lowest-lying 1 = resonance 
mass, assuming a resonance width of 400 MeV. 

Figure 12: Stability curves for determining the lowest-lying 1 = resonance 
mass, assuming a resonance width of 500 MeV. 

Figure 13: The s -dependence of the masses M a corresponding to (local) 
minima of the stability curves presented in Figs. 8-12. 

Figure 14: How improving the 4 square-pulse approximation affects M a : a 
comparison of three stability curves for the lowest-lying 1 = resonance 
mass that are obtained via (9) from weighting functions Wo^(M a , T, r) 
based upon truncation of the Riemann sum (11) at n = 4, 8, and 
14 pulses. A larger number of pulses corresponds to a more precise 
approximation of the Breit-Wigner shape than that presented in Fig. 
1. In all three curves, sq is assumed to be 1.63 GeV 2 , and the resonance 
width T is assumed to be 100 MeV. 

Figure 15: Stability curves for determining via (9) the lowest-lying 1=1 
resonance mass, denoted here as M s . The corresponding width of this 
resonance is assumed to be zero. As before, the stability curves are 
truncated to exclude the region in which the Borel mass scale M exceeds 



Figure 16: Stability curves for determining the lowest-lying 1 = 1 resonance 
mass, assuming a resonance width of 200 MeV. 
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Figure 17: Weighted least-squares fit of field theoretical content (circles) 
of Rq(t, s )/m 2 to the corresponding hadronic contribution (triangles) 
anticipated from the lowest-lying 1 = scalar resonance. Parameter 
values leading to this fit are m a = 0.95 GeV, r CT = 0, s = 3.3 GeV 2 , 
and g a /m 2 = 0.040 GeV 4 . 
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